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We report on observations of SN 1006 with Suzaku. We firmly detected K-shell emission from Fe, for 
the first time, and found that the Fe ionization state is quite low. The broad-band spectrum extracted 
from southeast of the remnant was well-fitted with a model consisting of three optically thin thermal 
non-equilibrium ionization plasmas and a power-law component. Two of the thermal models are highly 
overabundant in heavy elements, and hence are likely due to ejecta. These components have different 
' ionization parameters, riet ~ 1.4 x 10^° cm~^ s and n^t ~ 7.7 x 10® cm~^ s; it is the later one that produces 

Fe-K emission. This suggests that Fe has been heated by reverse shock more recently than the other 



^ I elements, consistent with a picture where the ejecta are stratified by composition with Fe in the interior. 

' On the other hand, the third thermal component is assumed to be solar abundance, and we associate 

c/3 , it with emission from the interstellar medium (ISM). The electron temperature (fcTg ~ 0.5 keV) is lower 

^ ' than that expected from the shock velocity, which suggests a lack of collisionless electron heating at the 
forward shock. The extremely low ionization parameter and extreme non-equilibrium state arc due to the 

04 ' low density of the ambient medium. 

\ Key words: Supernova Remnants: individual (SN 1006) — X-Rays: spectra 

■r^lj- ' 1. Introduction ct al. 1997; Winkler et al. 2005), and interpreted as being 

\ due to unshocked Fe in the interior of the remnant. These 

■ Based on the historical record, SN 1006 is widely re- results show that portions of the Fe-rich ejecta are still 

f — . [ garded to be one of the Galactic type la supernova rem- expanding freely, and have not yet been overrun by the 

. nants (SNRs) similar to Tycho's SNR (Schaefer 1996). reverse shock. However the amount of Fe inferred in these 

The current supernova (SN) models predict that Fe pro- studies (< 0.16 Mq) is much less than the amount (0.6- 

duction in type la SNe is far larger than that of core- 0.8 Mq) predicted to be produced in the thermonuclear 

' collapse SNe (e.g., Nomoto et al. 1984; Iwamoto et disintegration of a ~1.4 M© white dwarf. 



> 



JH \ al. 1999). Therefore, measuring the Fe abundance is one In the early phases of supernova remnant evolution, 

- - 1 of the essential clues to classify the SN type of any specific even the shock-heated plasma is far from thermal equi- 

remnant. The K-shell X-ray emission lines from ionized librium in terms of either the ionization or particle (elec- 

Fe offer the most direct information for abundance deter- tron and ion) temperatures. The ionization age, a key di- 

mination. In fact, strong Fe-Ka lines had been observed agnostic of the non-equilibrium ionization (NEI) state, is 

from Tycho's SNR (e.g., Hwang et al. 1998; Dccourchelle defined as riet, the product of the electron density and the 

et al. 2001). time since the gas was heated. Typically, riet is required 

In SN 1006, Vink et al. (2000) suggested the possible to be > 10^^ cm""^ s for full ionization equilibrium (Masai 

presence of an Fe-K emission line at 6.3±0.2 keV from 1984). The thermal X-ray spectrum of SN 1006 suggests 

the BeppoSAX MEGS spectrum. However, the detec- an ionization timescale of riet ~ 2 x 10^ cm^'^ s (Vink et 

tion was not clear, because the energy resolution of the al. 2000), which is far from the full ionization equilibrium, 

MEGS (a gas scintillation counter) was relatively poor and lower than nearly all other Galactic SNRs. The elec- 

(FWHM~500 eV @6 keV). Other recent X-ray missions, tron temperature at the northwest rim of the remnant has 

Ghandra and XMM-Ncwton, have not succeeded in de- been estimated to be kT^ ^ 0.7 keV (Long et al. 2003) or 

tecting an Fe-K line. Gold Fe in the interior of SN 1006 ~ 1.5 keV (Vink ct al. 2003). Since these values are much 

is known to exist based on ultraviolet absorption studies, lower than the ion temperature expected from the shock 

Blue and red-shifted Fe II absorption lines were detected in velocity of ~2980 km s~^ (Ghavamian et al. 2002), non- 

the spectra of background stars (Wu et al. 1993; Hamilton equilibration of ion and electron temperatures is present 
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Fig. 1. XIS-BI mosaic image of SN 1006 in the 0.55-0.6 keV 
band. Exposure and vignetting effects are corrected. The 
coordinates (R.A. and Dec) are referred to epoch J2000.0. 
The black square shows the FOV of the XIS on the SE region 
of SN 1006 where this paper is concentrated. The blue squares 
show the FOVs on the other three quadrants (NE, NW, and 
SW). 

in SN 1006 as well. 

The extreme non-equilibrium state and high shock ve- 
locity are likely due to the low density of the ambient gas, 
because SN 1006 {b = -1-14.6) is further from the Galactic 
plane than other Galactic SNRs, such as Gas A (6= —2.1) 
and Tycho's SNR (6 = +1.4). Although Kepler's SNR 
{b = 6.8) is nearly as far above the plane as SN1006 (when 
their respective distances are included) , Kepler appears to 
be evolving into a dense circumstellar medium. Therefore, 
the evolutionary state of SN 1006 may be the lowest of 
these youngest Galactic SNRs, even though the real age 
is ^3 times older than the others. 

In this paper, using the type la SNR, SN 1006, we con- 
sider the X-ray spectra to investigate this very early phase 
of SNR evolution. Suzaku (Mitsuda et al. 2007), with its 
high sensitivity for diffuse sources and spectroscopic res- 
olution in the ^^5-10 keV energy range, is well-suited for 
this study. 

We assume the distance to SN 1006 to be 6.8x10^^ cm, 
following Winkler et al. (2003). The errors quoted in this 
paper represent the 90% confidence level, unless otherwise 
stated. 

2. Observations and Data Reduction 

Four pointings were made of SN 1006 with Suzaku dur- 
ing the performance verification (PV) phase. Suzaku 
has one Hard X-ray Detector (HXD; Takahashi et 
al. 2007, Kokubun et al. 2007), and four X-ray Imaging 
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Fig. 2. Background-subtracted XIS spectra extracted from 
the whole SE quadrant (SN 1006 SE). The black and red 
points represent the FI and BI spectra, respectively. 

Spectrometers (XIS; Koyama et al. 2007) each placed in 
the focal plane of an X-Ray Telescope (XRT; Serlemitsos 
et al. 2007). This paper reports on the imaging and 
spectral results obtained with the XIS. The XIS con- 
sist of three Front-Illuminated (FI) GGDs and one Back- 
Illuminated (BI) GGD. The advantages of the former are 
high detection efficiency and low background level in the 
energy band above ~5 keV, while the latter has signifi- 
cantly superior sensitivity in the 0.2-1.0 keV band with 
moderate energy resolution. All four XRTs are co-aligned 
to image the same region of the sky. The XIS were op- 
erated in the normal full-frame clocking mode with the 
editing mode of 3 x 3 or 5 x 5. 

We employed the cleaned revision 0.7 data, and used 
the HEADAS software version 6.0.4 and XSPEG version 
11.3.2 for the data reduction and analysis. The X-ray 
data taken during low cut-off rigidity (< 6 GV) were ex- 
cluded so as to minimize any possible uncertainty in the 
non X-ray background (NXB). The total effective expo- 
sure times after these screening were ~50 ks, for each of 
the pointings. The response matrix files (RMF) and ancil- 
lary response files (ARF) were made using xisrmfgen and 
xissimarfgen (Ishisaki et al. 2007) version 2006-10-17. 

3. Overall Structure 

3.1. Thin- Thermal Plasma Map 

The field of view (FOV) of the XIS is - 18' x 18', and 
hence a four-point mosaic can completely cover SN 1006 
(about 30' diameter). Since we intended to study the ther- 
mal plasma, we searched for the region of the remnant 
with the most prominent thermal emission. We hence 
made a mosaic image in the OVII Ka line band (0.55- 
0.6 keV), because this line is a major component of the 
thermal plasma in SN 1006. The result is shown in fig- 
ure 1. We can see bright emission in this band in the 
northeast (NE) and southwest (SW) regions of the rem- 
nant. However, these are dominated by non-thermal emis- 
sion (Koyama et al. 1995), from which it is difficult to 
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Fig. 3. XIS intensity map at the Fe-Ka line (a: 6.33-6.53 koV band), from which the continuum flux at 6.1-6.3 keV band [shown 
in (b)] is subtracted. In both images, exposure and vignetting eff'ects are corrected. The data from the three FIs are combined. Two 
corners of the calibration sources are removed. The black squares indicate each FOVs of the XIS. The red ellipse shows the region 
where we extracted the spectra for a detailed analysis. 



extract thermal spectra. Other regions are dominated by 
thermal emission. In particular, as shown in figure 1, the 
southeast quadrant (here SN 1006 SE) is brighter than 
the northwest (NW). Therefore, our study concentrated 
on this region. The spectral results of non-thermal rims 
(NE and SW) are reported by Bamba et al. (2008), and 
those of the other regions (NW quadrant and the center 
of the remnant) will be reported in another paper in the 
future. The FOV of the XIS on SN 1006 SE is outlined 
by the black solid square in figure 1, which was observed 
on 2006 January 30 (Obs. 1D=500016010). 

3.2. Spectrum of the SN 1006 SE Region 

We extracted the spectrum of the entire SN 1006 SE re- 
gion, excluding the two corners in the FOV that contain 
calibration source emission. For the background, we used 
the North Ecliptic Pole (NEP) data (Obs. date = 2006 
February 10, Obs. ID = 500026010, Exp. time = -82 ks). 
Although we have been monitoring and correcting the in- 
crease of the Charge- Transfer Inefficiency (CTI) for the 
XIS (see Koyama et al. 2007), recovery of the energy res- 
olution cannot be made. However, we can ignore the dif- 
ference of the spectral resolution between the SN 1006 SE 
and NEP observations, because the latter observation was 
made only ten days after the former. To minimize the un- 
certainty due to the background subtraction, in particular 
that of NXB. we applied the same data-screening criteria 
to both the SN 1006 SE and NEP observations, and took 
background spectra from the same detector coordinates 
as the source regions after excluding point-like sources 
detected in the NEP data. The background-subtracted 
spectra are shown in figure 2. Since the data from the 



three FIs are nearly identical, we merged those individual 
spectrum to improve the photon statistics. 

As shown in figure 2, we found clear K-shell (Ka) lines 
from Ar. Ca, and Fc, for the first time. With a power-law 
plus Gaussian-line fit, we determined the line center en- 
ergy of the Fe-Ka to be ~6.43 keV. This energy constrains 
the Fe ionization state to be approximately Ne-like. 

3.3. Iron Line Map 

We show in figure 3a an image in a relatively narrow 
band (6.33-6.53 kcV) that contains the Fc-Ka line. This 
image was generated by subtracting the continuum fiux at 
energies of 6.1-6.3 kcV. (The image in this band is shown 
in figure 3b.) 

We can see that the Fe-Ka fiux is enhanced at the 
southern part of the remnant (outlined in red with a el- 
lipse), except for the NE and SW quadrants where the 
non-thermal emission is dominant. The mean surface 
brightness at 6.33-6.55 keV within the elliptical region 
is 8.5 (±0.5) X 10~^ photons cm~^ s~^ arcmin"^, while 
that outside it (only in the SE and NW quadrants) is 
4.6 (±0.3) X 10~^ photons cm~^ s~^ arcmin"^. In order 
to study the thin-thermal spectrum with the best S /N ra- 
tio for Fe-K line, we extracted the X-ray spectrum from 
within the elliptical region, excluding the corner of the cal- 
ibration sources. The background subtraction was made 
in the same way as that of the full-field spectrum. The re- 
sults are given in figure 8. Hereafter, all detailed analyses 
are made using this spectrum. 
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Table 1. The center energies and widths of the emission 
hues. 



Line 


Center energy* (eV) 


Widtht (eV) 


Mg-Ka 


1338 (1337-1340) 


< 5.4 


Si-Ka 


1815 (1813-1816) 


40 (38-42) 


S-Ka 


2361 (2355-2365) 


60 (54-65) 


Ar-Ka 


3010 (2991-3023) 


< 50 


Ca-Ka 


3692 (3668-3724) 


< 57 


Fe-Ka 


6430 (6409-6453) 


< 60 



* Errors (statistical only) are given in the parentheses 
(see text). 

^One standard deviation (Icr). 



3.4- Energy and Width of the Emission Lines 

In order to study the line features, we fitted the spectra 
extracted from the elliptical region with a phenomenolog- 
ical model; a power-law for the continuum and Gaussians 
for the emission lines. The best-fit central energies and 
widths for the emission lines are shown in table 1. Since 
the absolute energy calibration error is ±0.2% above 1 keV 
(Koyama et al. 2007), the relevant elements are uniquely 
identified from the best-fit central energy. These are also 
listed in table 1. Except for Mg, the central energies 
are significantly lower than that of the respective He-like 
Ka lines. We also note that the widths of the Gaussians 
identified with Si-Ka and S-Ka are significantly broader 
than the instrumental energy resolution at those energies. 
Figure 6a would help to compare the widths of the data 
with those of narrow lines. In subsection 4.3, we discuss 
this matter in detail. 

4. Spectral Structure in the Narrow Bands 

In order to study the plasma characteristics in relation 
to the various elements, we at first divided the spectrum 
into three representative energy bands: the 0.4-1.1 keV 
band for the O and Ne (light elements) lines, the 1.2- 
2.8 keV band for the Mg, Si, and S (medium elements) 
lines, and the 5-10 kcV band for the Fe (heavy element) 
line. With the fitting to these individual band spectra, 
we constrained the plasmas including light, medium, and 
heavy elements separately. 

4-.1. The 5-10 keV Band Spectra 

We fitted the 5-10 keV spectrum with a solar abun- 
dance (Anders and Grevese 1989) VNEI model. Since the 
NEIvers 2.0 plasma code does not include K-shell emission 
lines for ions below the He-like state, we reverted to the 
NEIvers 1.1 code. If we fix the Fe abundance to be solar, 
then the model cannot reproduce the Fe-Ka profile or flux, 
although the x^/d.o.f. of 25/28 is acceptable. Allowing 
the Fe abundance to be free results in a greatly improved 
best-fit xVd.o.f. of 10/27. The best-fit kT^, net, and Fe 
abundance are 3.8 (1.7-26) keV, 6.1 (2.2-10) xlO^ cm'^ s, 
and 6.5 (>1.9) solar, respectively. The Fe abundance is 
significantly enhanced relative to solar, even if the con- 
tinuum flux in this band is assumed to be composed of 
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Fig. 4. XIS spectra in the 5-10 kcV band with the VNEI 
modeh The black and red data points represent the FI and 
BI spectra, respectively. In the BI spectrum, the energy band 
above 8 keV is ignored. 

pure thermal emission. The best-fit spectrum is shown in 
figure 4. 

4.2. The 0.4-1.1 keV Band Spectra 

Figure 5 shows the 0.4-1.1 keV band spectra. This en- 
ergy band includes Ka lines of the light elements, with 
dominant emission from the oxygen K-shell line series. 
The spectra were fitted with a VNEI model, allowing the 
abundances of C, N, O, and Ne to be free parameters. 
Also, we let the abundances of Ca and Fe to be free, be- 
cause the L-shell lines of these elements fall into this en- 
ergy band (0.4-1.1 keV). Interstellar absorption was fixed 
to a hydrogen column density of 6.8 x 10^" cm~^ (Dubner 
et al. 2002). In the initial fits, we found a significant 
inconsistency between the FI and BI spectra near the en- 
ergy of the 0-edge (0.54 keV). This may be due to incom- 
plete calibration information for the contamination layer 
on the optical blocking filter (OBF) of the XIS (Koyama 
et al. 2007). Since the calibration of the OBF contamina- 
tion for the BI is more accurate than that for the FI, we 
decided to retain the BI data across this energy band and 
ignore the 0.5-0.63 keV band in the FI spectrum. The re- 
sults are shown in figure 5a, with the best-fit temperature 
and ionization parameter of kT^ = 0.58 (0.56-0.59) keV 
and net= 6.7 (6.6-6.8) xlO^ cm'^ s. 

In figure 5a, we find an apparent disagreement between 
the data and model in the '^0.7-0.85 keV energy band. 
Usually, this energy band is dominated by the iron L-shell 
line transitions of 2>s — 2p ('^730 eV) and 3d—2p ('^830 eV) 
from Fe XVII. Although the fluxes of these two lines in NEI 
models are nearly equal over a wide range of plasma condi- 
tions, extrapolation to the extremely low ionization states 
that we see here in SN 1006 is quite uncertain. In par- 
ticular, for the young Type la SNR E0509— 67.5, Warren 
and Hughes (2004) found that the Fe L-shell emission was 
dominated by an emission from a line feature near 0.73 
keV, which was not sufficiently strong in the spectral mod- 
els and had to be included as a separate Gaussian compo- 
nent. Based on the strength of other Fe-L emission lines 
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and the weakness of O emission in E0509— 67.5 the addi- 
tional hne feature at 0.73 keV was confidently associated 
with Fe-L emission. 

Although we cannot eliminate the possibility that the 
residual comes from faint Fe-L emission, here we con- 
sider a different explanation. The best-fit spectral model 
of the oxygen lines are shown in figure 5b, in which 
the black and red solid lines are K-shell transitions from 
O VII and O VIII, respectively. In most astrophysical plas- 
mas, the best-fit temperature of ~0.6 keV gives weak K- 
shcll lines of OVII compared to those of OVIII. On the 
other hand, in lower temperature plasmas (e.g., kT^ ~ 
0.1 — 0.2 keV) where the OviI Ka line is dominant, the 
line fiuxes decrease rapidly along the K-shell transition 
series (Ka,/3,7,i5,e,C, etc.). Therefore, K-shell lines in the 
higher transitions can be safely ignored, and hence conven- 
tional NEI codes do not include O VII K-shell transition 
lines higher than K(5. This is the reason that no oxygen 
line in the ~0.7-0.85 keV band is present in the best-fit 
model (figure 5b). In SN 1006, however, the ionization 
timescale is low, while kTe is moderate, so the fluxes of 
higher level K-shell transitions from OviI may be rela- 
tively strong and cannot be ignored. To account for them 
we added Gaussians at 714 eV, 723 eV, and 730 eV to rep- 
resent the OviI K5, Ke, and K^ lines, respectively. At a 
plasma temperature of 0.6 keV, the flux ratio of O VIII Lye 
to Ly5 predicted by the NEIvers 1.1 code is ~0.5 (see fig- 
ure 5b). In the absence of a detailed calculation we merely 
assumed that the flux ratio of the O VII lines follow the 
same pattern, namely Ke/K5 = K^/Ke = 0.5. The best- 
fit results with these additional lines, see figure 5c, are a 
reasonably good fit. The additional artificial lines are also 
shown as the dotted lines in figure 5b. The temperature 
and ionization parameter are almost the same as those ob- 
tained without the additional higher K-shell transitions of 
OVII. 

The best-fit intensity ratio of O VII K/3 to O VIII Lya 
is ~1.4, which is consistent with the result of the 
XMM-Newton/RGS observation of the X-ray knot lo- 
cafizcd at the NW rim of SN 1006 (Vink et al. 2003; 
OviI K/3/0 VIII Lya ~ 1.6). This agreement may justify 
our approach to include the K-lines manually. Also, this 
suggests that the O VII K-shell line series are key spectral 
components for all over the SN 1006 plasma. 

4.3. The 1. 2-2. 8 keV Band Spectra 

This band includes K-shell lines of medium-weight ele- 
ments (Mg, Si; and S). As shown in table 1, we see can 
significant broadening of the Si and S Ka lines. Therefore, 
the 1.2-2.8 keV band spectra cannot be reproduced with 
a single-component plasma model. In fact, we found sig- 
nificant disagreement between the data and model around 
the Si-Ka and S-Ka lines (figure 6a) in the one-component 
VNEI fits. The lack of data in the 1.83-1.85 keV band is 
due to current calibration errors of the XIS; there is a 
small gap in energy (~10 cV) at the Si K-edge, which is 
not implemented in the current response function. The 
apparent line broadening of the Si Ka line, however, is 
much larger (~40 eV) than this gap. Therefore, the dis- 
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Fig. 5. (a) XIS spectra in the 0.4-1.1 keV band fitted with a 
VNEI model. The black and red data points represent the FI 
and BI spectra, respectively, (b) Spectrum in photon space 
corresponding to the model used in (a), but only showing 
O lines. It is binned every 2 eV. The black and red lines 
represent O vil series and O VII series, respectively. The dot- 
tod lines arc additional transitions of K(5 (714 eV), Ke (723 
cV), and (730 cV) of Ovil added as separate Gaussians 
(see text), (c) Same spectra as (a), but for fits including 
Gaussians representing K5 (714 eV), Ke (723 eV), and K^ 
(730 cV) of O VII (see text). The residuals in the energy band 
0.7—0.85 keV seen in (a) are largely removed. 




Fig. 6. (a) XIS spectra in the 1.2-2.8 keV band fitted with the one-component VNEI model. The black and red represent the FI 
and BI, respectively, (b) Same to (a), but fitted with the two-component VNEI. 



agreement is not due to a calibration error, but is real. 

If the intrinsic broadening of the Si-Ka line (~ 40 eV; 
see table 1) is due to thermal Doppler broadening, the ion 
(silicon) temperature must be 13 MeV, which requires 
a shock velocity of ~ 1.5 x 10^ km s~^. Ghavamian et 
al. (2002) determined the shock velocity (from measure- 
ments of the Ha line width) to be ~2900 km s~^ in the 
NW portion of SN 1006. This is where the blast wave is 
interacting with significant amounts of interstellar matter, 
compared to other parts of the rim. In fact, the average 
size of the remnant (~15f5 radius) and the well-known 
age and distance allow us to determine a mean expansion 
speed for SN 1006 to be '-^lO'* km s^^. The true current 
blast wave speed in SN 1006 is almost surely bracketed by 
these two values. Furthermore, the reverse shock, which is 
likely to be the heating source for the Si that we see, typ- 
ically moves into the ejecta at only a fraction of the blast 
wave speed. Therefore, we consider it unlikely that the 
line broadening we see is due to thermal Doppler broad- 
ening. 

Another possible explanation proposes that the emis- 
sion consists of several thermal plasma components. We 
fitted the spectra with 2- VNEI models, in which we al- 
lowed for different kT^ values, but the same net between 
components. However, no combination of parameters 
would fit the data, even if we let the abundances of Mg, 
Si, and S be free parameters. Not only 2-VNEI, but 
also 2-VPSHOCK models failed to fit the data when the 
two components were forced to have the same ionization 
timescale. All of the models that we tried, but rejected. 



Table 2. Best- fit parameters and values of the spectral 
fitting in the 1.2-2.8 keV band with the various models. 



Model 


kTel 


kTe2 


net 


xVd.o.f. 




(keV) 


(keV) 


(cm~^ s) 




1-VNEI 


1.5 




4.5x10'-' 


848/351 


2-VNEI 


0.18 


20 


8.4x10^ 


706/346 


1-VPSHOCK 


1.6 




8.3x10^ 


726/351 


2-VPSHOCK 


0.36 


12 


1.0x10^" 


597/346 



are summarized in table 2. The addition of a power-law 
component to these models also did not help to improve 
the fittings. 

Next we tried two-component VNEI models with the 
same electron temperature, but different net values. The 
best-fit reduced x^/d.o.f. was greatly improved to 401/346 
(see table 2, for comparison). Models with different elec- 
tron temperature and ionization timescale {2-kTe and 2- 
net) gave no significant improvement of the reduced x^- 
We thus conclude that a 2-component VNEI model with 
different net values is necessary to fit the medium ele- 
ment plasma band. The best-fit parameters and spectra 
are given in table 3 and figure 6b, respectively. 

5. Model Fit of the Pull Band Spectra 

In section 4, we discuss how we separately derived spec- 
tral parameters for the three-energy bands that represent 
typical plasma conditions for the main light elements (O 
and Ne), medium elements (Mg, Si, and S), and heavy 
element (Fe). A summary is given in table 4. We used 
these parameters as the initial values to search for plasma 
parameters in the full energy band of 0.3-10 keV; the re- 
sulting overall best-fit parameters must be consistent with 
the data over the entire XIS energy band. 

Based on table 4, we assume that the model that de- 
scribes the full spectral range must consist of, at least, 
three plasma components: (1) high-fcTe (~1 keV) with 
high-TT-e^ (~ 10"'^° cm~^ s), (2) similarly high-fcTe with 
low-net (~ 10^ cm-3 s), and (3) low-fcTe (~ 0.6 keV) 
with medium-net (~ 7 x 10^ cm~'^ s). First, we consider 
whether or not the full band spectra can be well repro- 
duced with only these components. We also consider a 
broader range of kTe and net values for the Fe-K emis- 
sion. 

As shown in table 3, components (1) and (2) are sig- 
nificantly enhanced in Mg, Si, and S (relative to solar). 
In the present spectral fits, we made the assumption that 
these components are composed purely of metals without 
any admixture of hydrogen or helium. This assumption 



No. ] 

Table 3. 
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Best-fit parameters of the spectral fitting in the 1.2-2.8 keV band with 2-componcnt VNEI model with different riet values. 



Parameter 


Component 1 


Component 2 


iVn (cm-^) 


6.8x10^" 


(fixed) 


kTe (keV) 


1.1 (1.0-1.2) 


Mg 


4.2 (3.5-5.2) 


6.1 (3.5-13) 


Si 


19 (15-24) 


15 (13-24) 


S 


24 (19-34) 


23 (17-31) 


net (cm^'^ s) 


1.3 (1.1-1.7) xlO^'' 


7.9 (6.2-9.8) xlO^ 


nnUeV (cm^'^) 


2.8 (2.3-3.4) xlO^s 


1.9 (1.6-2.3) xlO^e 


xVd.o.f. 


401/346 


^ 1.16 



Table 4. Plasma components determined from the narrow band spectra. 



Band (keV) 


Major elements 


kTe (keV) 


net (cm s) 


Component* 


0.4-1.1 


0, Ne 


0.58 (0.56-0.59) 


6.7 (6.6-6.8) xlO'-* 


(3) 


1.2-2.8 


Mg, Si, S 


1.1 (1.0-1.2) 


1.3 (l.l-1.7)xl0i" 


(1) 






1.1 (1.0-1.2) 


7.9 (6.2-9.8) X 10^ 


(2) 


5.0-10 


Fe 


3.8 (1.7-26) 


6.1 (2.2-10)xl0^ 





*The plasma component identifications described in text. 



arose because hydrogen and helium emit only continuum 
emission in the X-ray band. Given the complexity of our 
spectral model, it is just not possible with the XIS data to 
obtain reliable estimates for the level of continuum emis- 
sion from the three plasma components separately. As a 
further assumption, we fixed the abundance of component 
(3) to be the solar value. (Anticipating the results given 
below, we found that most of the continuum emission 
above 1 keV actually comes from a hard power-law like 
spectral component.) Operationally, we fixed the oxygen 
abundance in both the pure-metal components (i.e., nos. 1 
and 2) to a large value (1 x lO''), and fitted for the abun- 
dances of other elements relative to oxygen. The C and Ni 
abundances were assumed to be the same as Ne and Fe, 
respectively. The emission measure (EM) in the plasma 
code is defined (in XSPEC) as n-aneV , even if the plasma 
is dominated by heavy elements. However, the oxygen 
density can be calculated as no = (8.51 x 10~'^rtH)xlO'^, 
where the numerical value is the solar abundance of oxy- 
gen from Anders and Grevesse (1989). Details of this 
method can be found in Vink et al. (1996). 

Since the NElvers 1.1 spectral model does not include 
K-shell emission lines from Ar and low ionization states of 
Ca, we added Gaussians at the energies of their Ka tran- 
sitions: 3.02 keV (Ar) and 3.69 keV (Ca) (see table 1), For 
component (1) we additionally freely varied the Ca abun- 
dance, which effectively allowed for the Ca L-shell lines to 
contribute at low energies (< 0.5 keV). For component (2), 
we tied the Ca abundance to that of Fe, which was allowed 
to vary freely. We also added Gaussians for the upper-level 
oxygen transitions at 0.714, 0.723, and 0.730 keV with the 
same reason as described in subsection 4.2. 

We fitted the FI data in the energy band from 0.33 to 
10 keV (0.33 to 8 keV for the BI data). We ignored the 
band below the C-edge at 0.28 keV, because of the large 
uncertainty in calibration due to contamination on the 
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Fig. 7. XIS spectra in the 0.33-10 keV band fitted with the 
three-component plasma model (see text) . The black and red 
points represent the FI and BI data, respectively. The model 
fails to describe the continuum emission above ~4 keV. 

OBF (see Koyama et al. 2007). Energies of 0.5-0.63 keV 
(only for FI) and 1.83-1.85 keV were also ignored for the 
same reasons as noted in subsections 4.2 and 4.3. Since 
the absolute gain of the XIS has an uncertainty of ±5 eV 
(Koyama et al. 2007), we allowed for a small offset in the 
photon-energy to the pulse-height gain relationship. 

With this model and assumptions, the best-fit x^/d.o.f. 
obtained was 1081/833. The fitted spectra are shown 
in figure 7. All line features are well fitted, even Fe-K, 
suggesting that including additional thermal components 
would not be justified. However, there is a systematic data 
excess in the continuum at energies ^ 4 keV. We include 
a power-law component to model this, since we cannot re- 
ject the possibility that there is a non-thermal component 
here that is similar to the bright NE and SW rims (albeit 
much less intense). The fit with the additional power-law 
was significantly improved with x^/d.o.f. = 996/831. The 
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Table 5. Best-fit parameters of the broadband spectral fitting. 



Component 


Parameter 






Value 


Absorption 


iVn (cm^^) 






6.8x10^*^ (fixed) 


VNEI 1 (Ejccta 1) 


kTe (kcV) 






1.2 (1.1-1.4) 




Abundance (10^ solar) 




C 


0.19* (0.08-0.30) 








N 


(fixed) 











1.0 (fixed) 








Ne 


0.19* (0.08-0.30) 








Mg 


4.1 (3.7-4.7) 








Si 


17 (15-20) 








S 


23 (20-27) 








Ca 


11 (9.6-13) 








Fe 


0.68^ (0.62-0.76) 








Ni 


0.68^ (0.62-0.76) 




Hgt (cm~^ s) 






1.4 (1.2-1.6) X 10^° 




nYiUfY (cm^'^) 






2.7 (2.4-3.0) X 10^1 


VNEI 2 (Ejccta 2) 


kTe (keV) 






1.9 (1.5-2.6) 




Abundance (10^ solar) 




C 


0.33* (0.23-0.45) 








N 


(fixed) 











1.0 (fixed) 








Ne 


0.33* (0.23-0.45) 








Mg 


3.0 (2.2-3.9) 








Si 


10 (8.6-12) 








S 


12 (9.1-16) 








Ca 


18§ (5.1-44) 








Fe 


18§ (5.1-44) 








Ni 


18§ (5.1-44) 




HfA. (cm~^ s) 






7.7 (6.7-9.2) X 10« 




n]ineV (cm""^) 






1.8 (1.3-2.5) X 10^2 


NEI (ISM) 


kTe (keV) 






0.51 (0.31-0.55) 




n^t (cm~^ s) 






5.8 (5.7-6.1) X 10^ 




niineV (cm""^) 






1.1 (1.0-1.2) X 10^*^ 


Power-law 


F 






2.9 (2.8-3.0) 




Norm" (photons cm^^ s 






6.3 (5.5-7.2) X 10-"^ 


Gaussian lines to complement the incomplete VNEI code (see text) 


Line 


Center Energy (keV) (fixe 


;d) 




Norm (photons cm ■^s ^) 


Vll-KS 


0.714 






6.4 (6.1-6.6) X lO-'^ 


VII-Ke 


0.723 






3.2 xlO""* * 


OVII-KC 


0.730 






1.6 xlO"'' ** 


Ar-K 


3.01 






5.7 (4.4-7.0) X 10-6 


Ca-K 


3.69 






2.4 (1.3-3.5) X 10-^ 


Gain 


Offset (eV) 






-4.8 


xVd.o.f. 








996/831 = 1.20 



*,*,"'' and ^ represent linked parameters. 

II The differential flux at 1 keV. 

#Fixed to 50% of the normalization of O VII-K(5. 

**Fixed to 50% of the normalization of O VII-Ke. 
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Fig. 8. XIS spectra in the 0.33-10 keV band which were extracted from the red elliptical region shown in figure 3. The black and 
red data points represent the FI and BI spectra, respectively. The solid lines, colored green, blue, light blue, orange, and magenta 
show the components of the ISM, Ejecta 1, Ejecta 2, Power-law, and additional Gaussians of the best-fit model for the BI spectrum. 
Note that the Gaussian like structures (magenta) at ~1.2 keV and '^1.9 keV are escape peaks of the added Gaussians at 3.02 keV 
and 3.69 keV due to the XIS response. 



best-fit parameters are given in table 5. As already noted, 
the gain may have an uncertainty of ±5 eV. The fitting 
process requires a gain offset of 4.8 eV for both FI and BI, 
within the allowable range. The best-fit model in the full 
energy band is shown in figure 8. 

6. Discussion 

In section 5, we analyzed the full energy band spec- 
trum extracted from the elliptical region shown in figure 3, 
and found that it could be described adequately with a 
parameterized model including three thermal plasmas in 
non-equilibrium ionization (VNEI 1, VNEI 2, and NEI, in 
table 5) and one power-law component (Power-law, in ta- 
ble 5). Two of the plasmas, VNEI 1 and VNEI 2, were as- 
sumed to have non-solar abundances, and the other, NEI, 
was assumed to have solar composition. In the following 
we discuss the implications of our spectral results, but 
the reader is cautioned that these results are subject to 
change if the assumptions we have made are changed. 

6.1. Origin of the Plasmas 

The NEI model, assumed to have solar abundances, pro- 
duces most of the observed low-energy X-rays, particularly 
the Ka lines from OVII, OVIII, and NelX (see figure 8). 
This component is rather uniformly extended over the en- 
tire remnant; as can be seen in the OviI line band map 
(figure 1). We associate this component with the swept-up 
ISM. 



The other plasmas, VNEI 1 and VNEI 2, with their non- 
solar elemental abundance ratios, are plausibly of ejecta 
origin. Since VNEI 1 has a larger ionization parameter 
than the other plasmas, we suggest that this plasma was 
heated by reverse shock in the early stage of remnant evo- 
lution (here Ejecta 1). On the other hand, VNEI 2 has 
an extremely low ionization parameter, and hence should 
have been heated much more recently (here Ejecta 2). We 
have firmly detected iron Ka emission for the first time. 
The low ionization state of Fe, plus its overabundance in 
the Ejecta 2 component, is generally consistent with the 
Type la origin of SN 1006. 

6.2. Relative Abundance in the Ejecta 

We compare our best-fit relative abundances with the 
predicted nucleosynthesis yield of the widely-used W7 SN 
la model (Nomoto et al. 1984). For Ejecta 1, although 
the abundances of C, Ne, Mg, Si, and S relative to O are 
broadly consistent with the W7 model, Ca and Fc fall far 
below their predicted values, as shown in figure 9a. In the 
case of Ejecta 2, on the other hand, the heavy elements, 
albeit with large errors, are broadly consistent with the 
abundance pattern from the W7 model, as shown in fig- 
ure 9b. These results, along with the difference in the ion- 
ization timescale between the components just mentioned, 
are consistent with a layed composition of the ejecta with 
the higher-Z elements more concentrated toward the cen- 
ter of SN 1006. 
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Fig. 9. Metal abundances as a function of atomic number derived from the spectral fitting. The data points of (a) and (b) represent 
those of Ejecta 1 and Ejecta 2, respectively. The solid lines show the abundance relative to oxygen calculated in the W7 model for 
a Type la supernova by Nomoto et al. (1984). 



For this simple interpretation, one expects the Fe line 
emission from the low-ionization timescale component to 
be spatially located interior to the Mg and Si emission 
from the high-ionization component, i.e., near the center 
of SN 1006. As shown in figure 3, though, the Fe flux 
appears to peak close to the southeastern rim. In sum- 
mary, our two-component spectral model for the ejecta 
is a highly simplified view of what is surely a complex 
multi-ngt and multi-fcTg structure that varies throughout 
the interior of the SNR. A longer observation, particularly 
of the inner regions of SN 1006, will help to improve the 
statistical accuracy of the Fe-K line flux, and allow us to 
study this apparent discrepancy in more detail. 

6.3. ISM Density and Ion Temperature 

The plasma parameters given in table 5 were derived 
from the southeast solid ellipse of figure 3, which is 
TT X 6fl X 4.'4 = 84 arcmin^ in area. This corresponds 
to 3.3 X 10"^* cm^, from which we estimated the emis- 
sion volume to be 1/ = (3.3 x 10^^)^/^ = 6.0 x lO" cm^. 
Therefore, the EM of the heated ISM derived from the 
spectral fits corresponds to an electron density of rig ~ 
0.15 /^°'^ cm^'^, where / is the filling factor. Since the 
age of SN 1006 is ^1000 yr, the ionization timescale can 
be roughly estimated as n^t ~ 4.7 x 10^ cm~^ s, which is 
almost consistent with the best-fit value of the ionization 
parameter of the NEI (ISM) component. Such a low ion- 
ization age suggests that the temperatures of the ions and 
electrons may also be far from equilibrium. According 
to equation (3) of Laming (2001), under the assumption 
that the initial ratio between kT^ and fcTu just at the 
shock front is the ratio of the electron and proton masses, 
the proton temperature is estimated to be, 

- ( s.Sxlotm-as ) ^ {oirtvT ^'^^^^ ' 

This result is consistent with that of the Ha observation of 
Te/Ta < 0.07 (Ghavamian et al. 2002). The extreme non- 
equilibrium state of the plasma is due to the low density 



of the ambient medium. According to equations (2) and 
(5) of Ferriere (2001), the density of Hi and Hll at the 
Z-height of SN 1006 {Z = 550 pc) is -0.03 cm'^. This 
is consistent with our estimate of the ambient density of 
nH/4'-0.03 cm-3. 

6.4- The Power-Law Component 

The low density of the ambient medium will allow the 
velocity of the shock front to remain high for a relatively 
long time. This may be one reason why this remnant 
shows such efficient particle acceleration up to very high 
energies (~ 10^^ eV), observed as power-law (synchrotron) 
emission from the NE and SW rims (Koyama et al. 1995). 
The photon index from these rims is F ~ 2.7-2.9. Now 
here, away from the bright rims, we also detect a power- 
law component in the hard energy band (the Power-law, in 
table 5) with a best-fit index of F ~2. 8-3.0, which is con- 
sistent with the NE and SW bright rim emission. Since 
the power-law emissions from the rims are most likely to 
be synchrotron X-rays, the power-law component newly 
found from the ellipse in figure 3 is also likely to be syn- 
chrotron X-ray emission. The surface brightness is esti- 
mated to be — 6 X 10~^^ ergs cm~^ s~^ arcmin"^ in the 
2~10 keV band. This value is ^ 50-times lower than the 
peak emission at the NE and SW rims. If similar to the 
NE rim, the power-law component should show thin fila- 
mentary structures (Bamba et al. 2003). However, XMM- 
Ncwton (Rothenflug et al. 2004) and Chandra (Hughes et 
al., in preparation) observations found little evidence for 
filamentary X-ray rims along the SE rim. 

7. Summary 

We have analyzed high-quality spectra obtained with 
Suzaku of a region in the southern part of SN 1006, se- 
lected because it is bright in Fe-Ka. The results and in- 
terpretations are summarized as follows: 

1. The spectrum can be described by a model with at 
least three NEI thermal plasmas and one power-law 
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component. 

2. The fits yield different ionization parameters of 
net 8 X 10^ cm"-^ s (low), 6 x 10^ cm~'^ s 
(medium), Ugt ^ 10^" cm~^ s (high). 

3. The low-Uet plasma are highly overabundant in 
heavy elements, in which we found Ka lines from 
Fe, for the first time. The abundance pattern is 
consistent with that of type la SN ejecta. 

4. The abundance of the medium-ngt plasma is as- 
sumed to be solar, and we associate this component 
with the shocked ISM. Although oxygen is not over- 
abundant, Ka lines of OVII and OVIII from this 
plasma appear to dominate the thermal emission 
from SN 1006. (However the complexity of our spec- 
tral model makes this claim less than definitive. We 
cannot exclude the possibility that the oxygen emis- 
sion from this ISM component actually comes from 
yet another ejecta component. Spatially resolved 
spectroscopy would help to resolve this concern.) 

5. The high-net plasma is overabundant in medium el- 
ements, like Mg, Si; and S, but heavy element like 
Fe. Our interpretation is that this plasma also has 
an ejecta origin, and the composition is dominated 
by lower atomic number species. 

6. Temperature equilibrium between ions and electrons 
has not yet been achieved. The proton temperature 
of the shocked ISM is ~15-times higher than the 
electron temperature. 

7. The spectrum from this region contains a power-law 
component. The photon index, F ~2.9, is similar to 
that of the northeast and southwest bright rims, sug- 
gesting that this component is also of synchrotron 
origin. 

The authors thank all member of the Suzaku team, es- 
pecially Y. Hyodo and H. Matsumoto. H.Y., S.K., H.N., 
and A.B. are supported by JSPS Research Fellowship for 
Young Scientists. 
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